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Abstract 

We discuss the stability of the Skyrmion solution in the presence of a finite isospin chemical 
potential \i. Solving numerically the mass of the Skyrmion as function of fx, we find a critical value 
fi c = 222.8 MeV where the Skyrmion mass vanishes. We compare the exact numerical treatment 
with an analytical discussion based on a special shape for the profile of the Skyrmion due to Atiyah 
and Manton. The extension of this ansatz for finite fi works quite well for /i < 121 MeV. Then, 
for small values of /i, where the analytical approach is valid, we consider the possibility of having 
an angular deformation for the Skyrmionic profile, which is possible for finite values of //. This is 
however, a small effect. Finally we introduce finite temperature corrections, which strength the 
instability induced by the chemical potential, finding the dependence of the critical temperature 
on \i. 
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I. INTRODUCTION 

The skyrmion picture has attracted the attention of many people as a possible way 
for understanding the hadronic dynamics, as well as the hadronic phase structure. The 
behavior of hadrons in the presence of a media, i.e. taking into account temperature and/or 
density effects, can be analyzed according to this perspective. For example, in a baryon- 
rich environment, it is found that coupling constants like g& and F v become quenched 
In reference P, I^J, which extends the construction of Atiyah and Manton in p, the 
authors discussed the stability of skyrmions under finite temperature conditions, showing 
the existence of a critical temperature T c . The skyrmion is not longer stable for T > T c . 

A similar kind of discussion has been carried out in the frame of the so called hybrid mod- 
els jf|. The construction of Hadrons as a core, described trough a chiral bag, surrounded by 
the tail of a skyrmion, has the same instability properties for a certain critical temperature. 
This can be interpreted as the occurrence of a deconfining phase transition. 

Recently Q, we have analyzed different topological structures in field theory when a 
finite isospin chemical potential (fi) is taken into account. In particular, we were able to find 
an upper bound fi c for the isospin chemical potential, such that the mass of the skyrmion 
vanishes. 

In this letter, we improve our previous discussion [7], showing that, in fact, the skyrmion 
develops an instability as function of fi, in the sense that the mass of the skyrmion, M = 
M(/i) diminishes and vanishes at a certain critical value fi c . We found fi c through an exact 
numerical analysis of the Skyrmion mass evolution, without referring to any particular radial 
profile. In a second stage, we compare the numerical results with an analytical treatment 
based on the Atiyah-Manton ansatz h\. We would like to emphasize that this profile is 
obtained from a construction based on SU(2) instantons in four dimensions, by computing 
Holonomies along lines parallel to the time axis. In this sense, this profile is definitely more 
fundamental than just an educated guess. The comparison between the numerical analysis, 
and the analytical procedure works quite well, in a wide range of fi up to a certain value of 
the order of fi ~ HOMeV, where in the analytical approach it is not longer possible to get 
a stable mass for the Skyrmion. 
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In addition, we propose an extension of the radial ansatz of the skyrmion when jti ^ 0, 
showing that a small angular deformation occurs for the minimum energy configuration, 
for small values of /i. This is done in the frame of a natural extension of Atiyah-Manton's 
ansatz. Finally we include temperature corrections, finding how the critical temperature 
depends on \x. 



II. ISOSPIN CHEMICAL POTENTIAL AND SKYRMION STABILITY 



The Skyrme lagrangian is 



+^Tr[{d»U)U\{d v U)tf}\ (1) 

where F n is the pion decay constant and e is a numerical parameter, following the nota- 
tion used in The introduction of an isospin chemical potential can be easily done by 
introducing a covariant derivative of the form P, Q| 
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Defining L v = (d u U)U\ the lagrangian density for finite \x becomes 



d v U -> D V U = d v U - i^[a\ U]g u0 - (2) 



£ M = ^Tr ^dJJdfrt + ^ [<T ~ Ua 3 tfa 3 ] j 

+ ^Tr^Luf-^ilQMuUL^Y (3) 

where g = a 3 — Ua 3 W and cr is the 2x2 unit matrix. 

With a little algebra, we can split the lagrangian into a zero and a finite chemical potential 
contributions 



fi=0 + 

/i 2 



64e 2 



Tr [g, L L 



(4) 
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£ M= o denotes the zero chemical potential term. 

The U field matrix can be parameterized in the standard way 



U = exp (— i£<? ■ h) = cos£ — i(a ■ n) sin£, (5) 

where a is the sigma matrix vector and h 2 = 1. This ansatz has a "Hedgehog" shape and 
the following boundary conditions have to be satisfied 



£(f) = f(r), n = r, 

e(o) = 7r, e(°o) = o. 



(6) 



The mass of the Skyrmion, for static solutions, develops a dependence on the Isospin 
Chemical potential as well as on the temperature. For the remaining of this section, we will 
concentrate on the T = scenario. If we define the non-dimensional variable f = eF n r, the 
mass of the Skyrmion will be given by 



Af„ = M„ =0 - 



-0—1 » r 

4e 3 F^ 2 32e 3 F- 4 ' 



where M At=0 is the zero chemical potential contribution 



M M=0 = 4vr / df 



Ae 



o 



2 ^1 W( « 



/■OO 

+47T / (if 

Jo 



sin 2 (6; 



X 



and ^2 and Z4 are the following non-dimensional integrals 



(7) 



(8) 



h 
h 



J rf 3 fTr [a - C/(7 3 C/V 3 ] , 



(9) 



We would like to emphasize that the chemical potential terms contribute with opposite 
sign. This implies that the solution becomes unstable above certain value of [i. 



Assuming a radial symmetric profile (£ = £( r ))> we can evaluate directly the traces in 
©, getting 



Tr [a - Ua 3 Wa 3 ] = 4 sin 2 9 sin 2 f, 



Tr [g, L t 



leading us to the following integrals 



47T 



32 sin 2 sin 2 f ( ^£ 



+ 4 



h — I dff 2 sin 2 £, 



32tt 



32 



sin 2 # sin £, 



sin 2 q-^l +-sin^ 
\ar / r z 



The variational equation for the profile according to (JJJ) is 



(10) 



(11) 



(12) 
(13) 



d 2 e , 



^ + 2sin^ - + ^ + sin2e f 



- 7 sm 2£ 

4 r l 



fi 2 r 2 sin 2 £ Z' ld 2 ^ I 



dA z + ld£ lsin2£ 2sin2£ 



0. 



(14) 



3 \ 2 (if 2 ' 2sin£ \df J f dr 4sin 2 £ f 

where fi = n/(eF n ). Equation f)14j) can be solved numerically for the profile for different 
values of fi (Figure^). Notice that the radio extension of the Skyrmion grows as a function 
of \x. 

In order to obtain the mass of the Skyrmion, the profile has to be inserted, numerically, 
in equation (|7Jl. Figure |21 shows the chemical potential dependence of the mass. The point 
where the mass vanishes corresponds to the critical value (i c = 222.8 MeV. This value does 
not depend on a model for the shape of £(r). It is a fundamental result associated to the 
Skyrmion picture. 

The same analysis can be done using a specific profile for £(r). Here we will use the 



ansatz introduced by Atiyah and Manton in 



which depends on a free parameter A: 



£ A (r) = Tr 



1 - 



Vr 2 + X 2 _ 



(15) 
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FIG. 1: Numerical solution for the profile £(r) for different values of \i. (// = 12.9 (MeV): solid 
line; fj, = 25.8 (MeV): dashed line; n = 38.7 (MeV): dotted line.) 
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FIG. 2: The numerical solution for the Skyrmion mass as funtion of /i. 

The variation with respect to A enables us to find the lower mass of the Skyrmion. This 
approach has proved to fit very accurately the numerical value of the mass of the Skyrmion 
in the /i = 0. 

The function £ can be scaled 



Mr) = £i(r/A), (16) 
d£x(r) 1 d£i(r/A) 

A d(r/X) ' 1 J 

in order to factorize A from the integrals in ()13|). Using y = r / X = f / X, with A = eF n X, we 
obtain 
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h = A 3 (yj J™dyy 2 vn 2 (Uy)/2) 
1.927 |'— ) A 3 



(18) 



Alf) [dyy*x 



sin 2 (6(l/)/2) (^) 2 + ^in 4 (6(l/)/2) 



= 3.728 I— A. (19) 

Defining the dimensionless quantities M M = M^e/F n and /i = fi/(eF w ), we find the 
Skyrmion mass as a function of the free parameter A 



M„ 



(4tt) 



0.504A + 4.254^ - 0.101u 2 A 3 
A 



- 0.311/1 2 A 



The A m which minimizes the expression takes the value 



(20) 



A 



0.523 — - 0.322 

jl 2 

1.038 



~4 



0.254- 8.511 /P + 0.097/1 



(21) 



where we have used the dimensionless parameter A and e = 5.45 according to |8j. In Figure 
Elwe show the dependence of A m on /i. In the limit fi — > 0, we recover the result A m = 2.904 
[5J]. Notice that we have an upper bound for jl = Jx c — > 0.173, where A m takes the value 
A m = 4.138. For bigger values of jl, X m becomes imaginary, implying, that the Skyrmion is 
no longer stable in the regime jl > jl c . The mass of the Skyrmion, as function of /i, shows a 
small decrease, about 8%. 

In Figure E] we show how the mass dependence on A evolves for three different values of 
the chemical potential. Note that for the jl = 0.926, the minimum disappears and we have 
an inflection point signalizing the end of the stable Skyrmion phase. 
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FIG. 3: The evolution of the \ m /F^ parameter as function of n/F n (The dot denotes the critical 
point). 
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FIG. 4: The Skyrmion mass as function of A for three values of fi (£i = 1.29(MeV): dot-dashed 
line; fi = 64.5 (MeV): dotted line; fi = 118.68 (MeV): continuous line ). 

In figure El we compare the numerical solution for the mass with our analytical result. We 
see that in a wide region, up to ji ~ 100 MeV, the discrepancy between both approaches is 
less than 2%, showing that the ansatz is a very good approximation in this region. 

Notice that the critical value /i c where the mass of the Skyrmion vanishes, is much lower 
than the one we obtained in our previous article Q|, using the Atiyah-Manton's ansatz with 
a fixed value for A. 
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FIG. 5: Comparison between the exact numerical solution (dashed line) and the Atiyah-Manton's 
ansatz (continous line) for the mass of the Skyrmion. 

III. ANISOTROPIC SKYRMION ANSATZ 

If fi 7^ 0, we may expect the occurrence of anisotropic stable skyrmions, i.e. solutions 
where the profile £ becomes also dependent on the azimuthal angle 9 (£ = £(r, 9)). The 
discussion that follows is based on an appropriate extension of Atiyah-Manton's construction 
in a region of /i, where we had almost a coincidence between the numerical and analytical 
solution. A possible simple ansatz which decouples the angular dependence in the integrals 
(J7J), is the following 

(22) 

where the function f{9) represents a small deformation with respect to the radial configu- 
ration. A natural candidate for f{9) is to take the first term in a Fourier expansion: 



6M) = 2tt 



1 - 



rf(9) 



f(6) = l + asm6 + bcos6. (23) 

We found that the cosine term does not contribute to dimmish the Skyrmion mass, so we 
can take 6 = 0. In this way, the energy equation is modified, due to the angular terms, 
according to 
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M^^(/ 1 + / 2 ) + i ^(/ 3 + « 

"*i.-=gU (24, 



4e 3 F n 32e 3 F ff 

In the previous expression, we have introduced the following set of integrals. 



h 

h 
h 

h 

h 
h 



2n I drr 2 d9 sin9 
2n I drr 2 d9 sin 9 



(Mi) 2 + ^{doti) 



COS^i 



it / drr d9 sin 9 



1 — COS^i 



x {diixdjh - djiidih) 



2n J drr d6 sin 9 (1 — cos£i) (difi x dj-hy 

8tt J drr 2 d9 sin 3 9 9 sin 2 (fj/2) 

2tt / drr 2 dfl sin fl{8 sin 2 flsin 2 (£i/2) 



x 



(<9,6) 2 + ^(^i) 



32 • 2. • 

H — =- sin 9 sm 



l (a/2)}. 



(25) 



In Fig |U1 we have plotted the dependence of the deformation parameter a(/i), showing 
that there is a small shift in the shape of the Skyrmions. The criteria for determining the 
dependence of a(/x) is to minimize equation Such deformation is possible due to the 

special character of the Skyrmionic solution which is a topological object that entangles the 
isospin with the spatial degrees of freedom and therefore, a finite isospin chemical potential 
has a nontrivial consequences. 
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FIG. 6: The dependence of the deformation parameter a as function of \i. 
IV. FINITE TEMPERATURE EFFECTS 

Finally, we will include temperature effects in our discussion, following closely [3], using 
their thermal profile for the Skyrmion 



ff(r) = 2n 



1 - 



r + |A 2 (k coth(fix) — 1/r) 
■ 2 + tK 2 A 4 + kt\ 2 coth(fi;r) 



(26) 



where k = 2nT. In the limit k — >• 0, we recover equation f)15|) . Such profile was derived 
through an extension to finite temperature of the instanton construction by Atiyah and 
Manton 
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FIG. 7: The dependence of the critical temperature on \i. 
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In this part of the analysis we will not take into account possible angular deformations of 
the Skyrmion, since temperature does not affect the shape of the Skyrmion. Also the onset 
of instability due to chemical potential effects will occur for lower values of \i and therefore 
we can disregard it. 

The procedure we follow is to insert the thermal Skyrmion ansatz (|26|) in equation (J7|). 
Then we minimize with respect to the parameter A, getting for a certain /z, the critical 
temperature T c where the solution cannot attain a minimum, becoming unstable. 

In figure (J7J) we show how T c depends un /i. It is clear that T c diminishes with /z. 
Both parameters, temperature and chemical potential, cooperate in order to increase the 
occurrence of instabilities in the Skyrmion solution. 

In this letter we have discussed the stability of the Skyrmion configuration in the presence 
of finite isospin chemical potential (/z), showing the existence of a critical value /i c where 
the mass vanishes. Since the idea behind the Skyrmion approach is to have an effective 
representation of baryons, the existence of such phase transition is quite interesting. There 
are in the literature other discussions of phase transitions induced by density and tempera- 
ture. For example, in the frame of chiral Lagrangians, a detailed analysis of the phases of 
pion systems |ll[ has been carried out, showing the existence of phase transitions with the 
same kind of behavior for the masses that we have presented here, in a complete different 
approach. 
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